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KNO scaling, i.e. the collapse of multiplicity distributions P„ onto a universal scaling curve manifests when P„ 
is expressed as the distribution of the standardized multiplicity (n — c)/A with c and A being location and scale 
parameters governed by leading particle effects and the growth of average multiplicity. At very high energies, 
strong violation of KNO scaling behavior is observed (pp) and expected to occur (e"'"e~). This challenges one to 
introduce novel, physically well motivated and preferably simple scaling rules obeyed by high-energy data. One 
possibility what I find useful and which satisfies the above requirements is the repetition of the original scaling 
prescription (shifting and rescaling) in Mellin space, that is, for the multiplicity moments' rank. This scaling 
principle will be discussed here, illustrating its capabilities both on model predictions and on real data. 
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1. THE EARLY YEARS 

1.1. Polyakov and Koba-Nielsen-Olesen 

In the Millennium Year we celebrated the 30th 
anniversary of the famous Eq. (1), the basic re- 
sult of Polyakov and of Koba, Nielsen and Olesen 
concerning the asymptotic behavior of the multi- 
plicity distributions ■ These authors put for- 
ward the hypothesis that at very high energies s 
the probability distributions Pn{s) of producing 
n particles in a certain collision process should 
exhibit the scaling (homogeneity) relation 

{n{s)) \{n{s))J 

as s — > oo with (n(s)) being the average multi- 
plicity of secondaries at collision energy s. This 
so-called KNO scaling hypothesis asserts that if 
we rescale P„(s) measured at different energies 
via stretching (shrinking) the vertical (horizon- 
tal) axes by (n(s)), these rescaled curves will co- 
incide with each other. That is, the multiplic- 
ity distributions become simple rescaled copies 
of the universal function '4'{z) depending only on 
the scaled multiphcity z = n/{n{s)). In the pic- 
turesque terminology of Stanley Q the rescaled 



data points P„ (s) measured at different energies s 
collapse onto the unique scaling curve ij{z). After 
Koba let me illustrate this schematically: 
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The data collapsing behavior is a fundamental 
property of homogeneous functions. According 
to Eq. (1), at very high energies the multiplicity 
distributions Pn(s) are homogeneous functions of 
degree — 1 of n and (n(s)). Homogeneity rules 
play a central role in the theory of critical phe- 
nomena m . Near the critical point of a physical 
system the thermodynamic functions exhibit ho- 
mogeneous form which implies the existence of a 
law of corresponding states: using a suitably cho- 
sen scaling transformation it is possible to bring 
different states of the same system to coincidence 
and thus to compress many experimental or the- 
oretical results into a compact form. The KNO 
scaling hypothesis Eq. (1) is just such a law of 
corresponding states. 

In the celebrated paper |^ Koba, Nielsen and 
Olesen showed that the validity of Feynman scal- 
ing is a sufficient condition for Eq. (1) to hold. 
It is worth mentioning however that the mathe- 
matical rigor of their derivation was criticized by 
several authors. Moreover, the precise finding of 
Ref. H] was the energy independence of the mo- 
ments Cq = {n'^{s)) / {n{s))'^ in the limit s ^ oo. 
The scaling hypothesis Eq. (1) was formulated 
with the additional assumption that the moments 
Cq = /p z'^ip{z) dz determine uniquely the scal- 
ing function il}{z) . Polyakov arrived at the asymp- 
totic multiplicity scaling law Eq. (1) two years 
earlier than Koba, Nielsen and Olesen by formu- 
lating a similarity hypothesis for strong interac- 
tions in e^^e^ — > hadrons annihilation; see later. 
Somehow this fundamental and elegant work was 
almost completely ignored in the early 70s and 
Eq. (1) became known as the KNO scaling hy- 
pothesis in the high-energy physics community. 

1.2. Modifications 

The scaling relation Eq. (1) can hold only ap- 
proximately since multiplicity n is a discrete ran- 
dom variable whose stretching or shrinking by a 
scale factor leaves the probabilities P„ unaltered. 
The proper meaning of Eq. (1) is that with in- 
creasing collision energy s the discrete multiplic- 
ity distributions P„ can be approximated with 
increasing accuracy by a continuous probability 
density function f{x) via Pn ^ f{x — n) (KNO 
prescription) or by Pn ~ J^Zn^^ fi^) '^^ (called 



KNO-G "scaling") where, in the most popular 
modification of the original scaling rule, f{x) has 
the generic Czyzewski-Rybicki form [^,^ 

with scale parameter A > and location parame- 
ter c. Data collapsing behavior is observed if the 
only s-dependent parameters of the approximate 
shape function f{x) are c and A, i.e. if f{x) de- 
pends on collision energy only through a change 
of location and scale in x. The variation of A with 
increasing s reflects the growth of average multi- 
plicity, whereas the s-dependence of c is usually 
associated with leading particle effects (so-called 
AKNO or KNO-a scaling). Other physical expla- 
nations arc also possible, see e.g. the contribution 
of Marek Ploszajczak on A-scaling. 

2. NEW SCALING PRINCIPLE 

Taking into account a possible location change 
(shift) in multiplicity rarely proved to be suffi- 
cient to restore the data collapsing behavior of 
Pn when the original scaling law Eq. (1) is vio- 
lated at very high energies (hence we use c — 0). 
Nevertheless, the principle of scale and location 
change in the shape of multiplicity distributions 
can be extremely powerful through the following 
generalization of the KNO scaling rule. Consider 
that the multiplicity distributions, more precisely 
the continuous probability densities f{x) approx- 
imating Pn generate s-dependent shifting and 
rescaling in Mellin space. The Mellin trans- 
form of a probability density f{x) is defined by 
M{f{x);q} ~ x''~^ f{x) dx and it provides 
the (g — l)th moment {x''~^). Via the functional 
relation 




one can introduce translation and dilatation in 
the moments' rank q by performing the transfor- 
mation f{x) — > x^ f[x^^^ of the probability den- 
sities approximating the shape of P„. In the fol- 
lowings I shall illustrate the utility of performing 
a location or scale change in Mellin space; for fur- 
ther details, see Refs. 01 . 
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3. TRANSLATION IN M-SPACE 

3.1. Motivation 

There are several possible dynamical explana- 
tions of a shift in Mellin space. The basic idea 
is simply the following. Recall that in Eq. (1) 
the normalization of the KNO function reads 
ip{z)dz — zip{z)dz = 1, that is, (z) = 1. 
Assume that violation of KNO scaling is observed 
and we measure s-dependent "scaling" functions 
tp{z). In order to arrive at data collapsing be- 
havior of Pn the simplest possibility is to rescale 
the function ip{z) = zip{z) by its first moment 
given by C2. In terms of P„ and the moments 
(n'^) this corresponds to rescaling the distribu- 
tion Pn,i = nPn/ (n) by its mean {n) 1 = (n^) / (n) . 
The distribution Pn,i is known in the mathemat- 
ical literature as the first moment distribution of 
Pn. If the rescaling of Pn,i by (n)i is not suf- 
ficient to arrive at data collapsing behavior, one 
can try the same recipe for the higher-order mo- 
ment distributions Pn^r = Pn/ (n^)- The gen- 
eralization is straightforward. Let me show you 
now an example when Eq. (1) is clearly violated, 
but performing a shift in Mellin space can restore 
the collapse of different P„ curves onto a single 
scaling curve. 

3.2. Intermittency and multifractality 

The most obvious example is intermittency. As 
it turned out during the past 15 years, multipar- 
ticle production often yields scale-invariant den- 
sity fluctuations ||l|-|l|l. This can be observed 
through the power-law scaling Cq cx d^'^i of the 
normalized moments Cq = (n')/(n)'^ of Pn{S) as 
the bin-size 6 in phase-space is varied (for clarity, 
we neglect here the influence of low count rates). 
Assuming that the fluctuations show monofrac- 
tal structure, the so-called intermittency expo- 
nents ipq are given by (/s^ = (/92('Z ~ 1) and the 
anomalous fractal dimensions Dq — ipq/ (q — 1) 
are g- independent, Dq = D2. The normalized 
moments Cq of Pn{5) take the form 

Cq^Aq[C2]'>-^ for q>2 (4) 

with coefficients Aq independent of bin-size 5 [|l2| . 
Of course Eq. (1), with s ^ 5, is violated since 
we observe C2 c>c S^'f^ . 



In the restoration of data collapsing behavior 
of Pn for self-similar fluctuations the basic trick 
is the investigation of the higher-order moment 
distributions defined before. Their moments are 
{'nfl)r — {n'^^'^) / {rf) , that is, the moments of the 
original Pn are transformed out up to rth or- 
der via performing a shift in Mellin space, see 
Eq. (3). For r = 1 the normalized moments of 
the first moment distribution Pn^i are found to 
be Cq^i = Cqj^i/[C2\ in terms of the original Cq 
and comparison to Eq. (4) yields Cq^i — Aq+i 
for monofractal multiplicity fluctuations. Since 
the coefficients Aq are independent of bin-size 5, 
we see that monofractality yields not only power- 
law scaling of the normalized moments of P„ 
but also data collapsing behavior of the first mo- 
ment distributions Pn,i measured at different res- 
olution scales 5. Increasing the (possibly frac- 
tional) rank r of the moment distributions allows 
the restoration of data collapsing behavior in the 
presence of increasing degree of multifractality of 
scale-invariant multiplicity fluctuations [0,D . The 
effect of low multiplicities (Poisson noise) can be 
taken into account via the study of factorial mo- 
ment distributions Pn,r = n^^^ Pn / (n^^^) and their 
factorial moments. 

3.3. Ginzburg-Landau phase transition 

In recent years considerable interest has been 
devoted to the Ginzburg-Landau theory of the 
phase transition from quark-gluon plasma to 
hadronic matter ||l^,|l^ . Instead of a strict power- 
law scaling Fq cx S^''^'' of the normalized fac- 
torial moments, a so-called P-scaling behavior 
Fq cx [F2f-' is obtained with Pq = {q - . For 
second-order transition ly ~ 1.304 and in case of 
a first-order transition 1/ — 1.45. The F-scaling 
rule of the above form makes difficult to iden- 
tify the family of probability laws P„ belonging 
to the model. But considering P„_i instead of P„ 
the corresponding factorial moments turn out to 
be Fq^i (X [P2]''' with Pq ^ q" - q and hence 

Fq,l^[F2^iP With = 

This form of P-scaling is the familiar log-Levy 
law |]l^ with stability index < ly <2. Thus, via 
moment shifting, we succeeded in identifying P„^i 
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for the Ginzburg-Landau formalism with a strict 
bound on the essential parameter v of the model. 
The theoretically relevant values are within the 
allowed range. Although this particular example 
is not related directly to data collapsing behav- 
ior of multiplicity distributions, it illustrates the 
utility of the moment distributions P„,r oc rfPn 
which correspond to translation in Mellin space. 

4. RESCALING IN M-SPACE 

4.1. Motivation 

The most important reason of a possible change 
of scale in Mellin space comes from QCD. In 
higher-order pQCD calculations, allowing a more 
precise account of energy conservation in the 
course of multiple parton splittings, the natu- 
ral variable of the multiplicity moments is the 
rescaled rank q'f instead of rank q itself ||20|-[2^ 
with 7(0;^) being the multiplicity anomalous di- 
mension of QCD. Because of the running of the 
strong coupling constant as, it is inevitable to 
adjust an energy dependent scale factor in Mellin 
space if one wants to arrive at data collapsing 
behavior of the multiplicity distributions P„(s). 

4.2. Running coupling 

In the very first paper predicting the multiplic- 
ity scaling law Eq. (1), Polyakov constructed a 
model for the short-distance behavior of strong 
interactions based on the hypothesis of asymp- 
totic scale- and conformal invariance |jl|. When 
applied to the process of e+e" annihilation to 
hadrons at asymptotically high collision energies, 
the model predicts that multihadron production 
goes in a cascade way. First, a few heavy vir- 
tual objects, called jets, are formed. These are 
repeatedly diminished in a branching process un- 
til the masses of the subsequent generation jets 
become at some very late stage comparable to 
the hadronic masses. This cascading mechanism 
gives rise to multiplicity distributions P„(s) sat- 
isfying Eq. (1) and the scaling function ifi{z) be- 
haves as ip{z) cx a{z) exp(— z^) with /i > 1 and 
a{z) being a monomial. Up to a scale factor, 'ip{z) 
is a negative binomial type scaling function, that 
is, a gamma distribution ||l^ in the rescaled and 
power-transformed multiplicity z'^. 



It is worth emphasizing that the above result 
was achieved years earlier than the emergence 
of quantum chromodynamics. Interestingly, the 
predictions of QCD obey some important simi- 
larities to the findings of Polyakov: 

1 ) Hadronic jets produced in hard processes are 
expected in QCD to be self-similar objects com- 
posed of subsequent parton branching decays . 
A highly virtual quark or gluon decays into secon- 
daries which are less off-shell and less energetic. 
Each of these intermediate decay products splits 
into novel ones. The subsequent decay chain steps 
form a scale-invariant branching process which 
continues until the virtual mass of the quanta be- 
comes approximately 1 GeV. 

2) The fragmentation of partons in QCD leads 
to the typical characteristics of branching pro- 
cesses, in particular, to long-range correlations 
and to the validity of KNO scaling p3|. Further, 



it was shown 1 19| that KNO scaling in e+e anni- 



hilation is the consequence of the scale-invariance 
(asymptotic freedom) in the presence of gluon 
self-coupling. It is thus the manifestation of the 
essence of the non-abelian gauge theory of strong 
interactions. 

3) Taking into account the higher pQCD ef- 
fects responsible for energy-momentum conserva- 
tion in parton jets, the KNO scaling function in 
e^e^ annihilation was shown ||20|| to be identical 
to Polyakov's form: 

= AAz^'=-i exp ( - [Dzf) (5) 

where k = 3/2, M = i^D'''' /T{k), fi = (1-7)-^ « 
5/3 and D is a scale parameter depending on 
7, with 7 being the QCD anomalous dimension, 
7 « 0.4 at LEPl. Hence, the conservation laws 
strongly reduce the width of the KNO scaling 
function as compared to the exponential fall-off 
obtained by lower-order pQCD calculations. The 
experimental data at ^/s — 91.2 GeV (dots) con- 
firm the prediction (curve) as illustrated in the 
figure a) below (with the exception that k turns 
out to be larger, fc « 5). 

The crucial difference between the QCD pre- 
diction and Polyakov's result lies in the fact that 
the scaling function ip{z) given by the Polyakov- 
Dokshitzer form Eq. (5) depends on collision en- 
ergy s in QCD. That is, QCD predicts violation of 
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KNO scaling in e+e^ annihilation which becomes 
clearly visible at higher energies. This effect is 
due to the running of the strong coupling constant 
Us, hence 7 ^ as s ^ oo, i.e., /Lt ^ 1 asymptot- 
ically ll^. Therefore the tail of tpi^) widens with 
increasing s which gives rise to the KNO scaling 
violation pattern shown in the figure b) below. 
Asymptotically the exponential fall-off predicted 
by the double logarithmic approximation is recov- 
ered (solid curve) i.e. the negative binomial type 
scaling form of P„(s). 



0.01 - 



0.0001 




0.01 - 



0.0001 




1 2 
z 



=1 



0.01 



0.0001 




1 



Thus we have learned that the pre-QCD and 
QCD-based descriptions of e+e^ annihilation to 
hadrons make essential use of self-similar branch- 
ing processes, further, both predict the rescaled 
and power-transformed multiplicity variable z'^ 
to be gamma-distributed. But in the pre-QCD 
approach the exponent fj, is independent of col- 
lision energy s and KNO scaling holds valid, as 
expected for scale-invariant branching with fixed 
coupling. On the contrary, the QCD-based calcu- 
lations predict violation of KNO scaling since fj, 
varies with collision energy s, due to the running 
of as- Note however that in the latter case data 
collapsing can be restored in a simple way using 
logarithmic scaling variable. For the Polyakov- 
Dokshitzer form of ipiz) given by (5) we have 

ipix) = ^exp (fc^x - e'"")/r(fc), X = ln{Dz). 

Since only the exponent fi and scale parameter 
D of Eq. (5) are expected to depend on collision 
energy s through the variation of 7 (as), KNO 
scaling violated by QCD effects is recovered by 
plotting /i~^'(/'(/ix) as displayed in the figure c). 
The scale change in logarithmic multiplicity is 
governed by the QCD multiplicity anomalous di- 
mension. This particular type of data collapsing 
of the multiplicity distributions -Pn(s) is called 
log-KNO scaling since we have the behavior 
of type Eq. (2) except that distributions of the 
logarithmic multiplicity satisfy it. 

4.3. Multiplicative cascades 

It may turn out that the log-KNO scaling law 
has ubiquitous appearance in multiparticle dy- 
namics. Random multiplicative cascades play 
a distinguished role in the dynamics underlying 
multihadron production, both in soft and hard 
processes. The multiplicative cascade models like 
for example the a-model and p-model are in the 
focus of interest since the pioneering work of 
Bialas and Peschanski , see also 1 12 1 for a re- 
view. Recently, Frisch and Sornette developed a 
theory of extreme deviations |2^ which predicts 
for multiplication of random variables the generic 
presence of stretched exponential distributions 
f{x) cx exp(— Ax^) with ^ < 1. To be specific, 
consider the product Xjv = mim2 • • • mjv of in- 
dependent random variables distributed iden- 
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tically according to the probability density p{mi). 
The above product can arise in a cascade process 
if primary entities (e.g. particle density) of initial 
size So are repeatedly diminished in random pro- 
portions, so that after A'^ cascade steps we have 
the product sn = ■ • ■ niN- 




In ref. | |23[ it is shown that under mild regularity 
conditions the distribution of the product variable 
Xjsi exhibits the asymptotic behavior 

Vn{X)^[p{X^''')]^ for X^oo (6) 

and for finite N . Frisch and Sornette call atten- 
tion to the intuitive interpretation of (6): the tail 
of Vn{X) is controlled by the realizations where 
all terms in the product are of the same order 
hence Vn{X) is, to leading order, just the prod- 
uct of the N densities p, each having the com- 
mon argument X^^'^ . When p{x) is chosen to 
be oc exp(— Ax") with a > 0, then Eq. (6) leads 
to stretched exponential tails cx exp(— AiVx"/^) 
for large A^, with stretching exponent a/N < 1. 
Expecting the cascade depth N to be an increas- 
ing function of collision energy s, Eq. (1) breaks 
down but data collapsing can be recovered in the 
log-KNO manner discussed previously. 

5. SCALING AT TEV ENERGIES 

The most exciting and challenging task in de- 
veloping novel scaling relations is testing them 
on real data. In case of multiplicity distributions 
the KNO scaling laws Eqs. (1-2) are known to be 
strongly violated above ISR energies. This viola- 
tion is most visible for the multiplicity data mea- 
sured by the E735 Collaboration (H) see also 
for more details. In the foUowings I provide a 
very brief summary of some preliminary results 
of a study of the E735 data. 



5.1. Analysis of E735 multiplicities 

The E735 Collaboration recently published the 
full phase-space multiplicity distributions in pp 
and pp collisions at cm. energies = 300, 546, 
1000 and 1800 GeV at Tevatron ||. The Pn 
curves corresponding to -y/s = 300 and 1800 GeV 
are displayed in the following two figures. 




50 100 150 200 250 n 

Pp n I I I I I 




50 100 150 200 250 n 



It is apparent that bimodal shape of the P„ 
curves arises at TeV energies, obeying a (not too 
pronounced) shoulder structure like at SPS. It is 
argued |^,^ that the low-multiplicity regimes 
are affected mainly by single parton collisions and 
exhibit KNO scaling, whereas the large-n tail of 
the distributions is influenced more heavily by 
double parton interactions and violate KNO scal- 
ing substantially. Another possible explanation of 
the observed P„ shape is the weighted superposi- 
tion of a soft and semihard component, the lat- 
ter one being dominated by minijet production. 
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These components are responsible of the KNO 
scahng and KNO-violating regimes of P„ ■ 
The sohd curves in the previous two plots rep- 
resent the Polyakov-Dokshitzer parametrization 
Eq. (5) fitted to the KNO-violating z > 1 com- 
ponent of the E735 data. The quality of fits is 
reasonably good (for all data sets) as can be de- 
picted from the figures. The fits correspond to 
k = \ kept fixed and /i decreasing from /i ^ 2.2 to 
fj, ~ 1.7 as one goes from ^/s = 300 GeV to 1800 
GeV. Our findings indicate that the KNO violat- 
ing component of the E735 multiplicity distribu- 
tions can be collapsed onto a single scaling curve 
in the log-KNO manner discussed previously. It 
is illustrated below. 
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In the upper figure the ^/s = 300 and 1800 GeV 
data sets are displayed only, whereas in the lower 
plot all the four data sets measured by the E735 
experiment. The collapse of the data points onto 
a unique scaling curve is apparent. This very pre- 



liminary investigation suggests that the pattern 
of KNO scaling violation by double parton col- 
lisions, or by semihard processes dominated by 
minijet production, is similar to that of e^e~ an- 
nihilations where log-KNO behavior is expected 
to arise as well at high energies (see before). But 
note that in parameter k there is an order of mag- 
nitude difference between e^e~ and pp collisions, 
that is, the scaling functions have different shapes 
in the two reactions. 

One of the earliest reviews of Polyakov's pi- 
oneering work concluded with the following spec- 
ulative note: "Could it be that as in the the- 
ory of critical phenomena, while the basic physics 
(Hamiltonians) between e"*"e~ annihilation and 
hadronic multiparticle production could be en- 
tirely different, multiplicity scaling may yet still 
be a universal feature shared by both?" In the 
light of our findings the answer to this question 
is yes but the scaling behavior foreseen 30 years 
ago arises in the logarithm of multiplicity instead 
of multiplicity itself. 

6. SUMMARY 

In the past 30 years, data collapsing of multi- 
plicity distributions P„ according to the generic 
KNO scaling relation Eq. (2) was one of the most 
extensively studied topics in soft physics. Unfor- 
tunately, more and more experimental and the- 
oretical results confirm that Eq. (2) is too sim- 
ple to be true. The information content of P„ 
can be represented in an equally convenient man- 
ner making use of the multiplicity moments. For 
example, validity of Eq. (2) corresponds to con- 
stancy of the normalized central moments for ap- 
propriately chosen c and A. What can we do if 
the KNO scaling rule (2) breaks down? How data 
collapsing behavior of P„ can be restored, if it 
is possible at all? The approach presented here 
is extremely simple: let's perform the original 
scaling transformation (shifting and rescaling) in 
the multiplicity moments' rank, that is, in Mellin 
space, rather than in multiplicity. 

The functional relation Eq. (3) tells how things 
transform at the level of P„. Translation in 
Mellin space (moment shifting) corresponds to 
changing from Pn to the moment distributions 



8 



Pn,r oc n^Pn- This makes possible to arrive at 
data collapsing for intermittency phenomena: e.g. 
Pn,i scales for monofractal fluctuations, P„.2 for 
bifractals and so on. Dilatation in Mellin space 
(rescaling of moments' rank) corresponds to the 
transformation Pn Pm with m = n^^. There- 
fore data collapsing arises in rescaled logarith- 
mic multiplicity (log-KNO). Such behavior is ex- 
pected to occur for phenomena governed by run- 
ning coupling effects or in multiplicative cascades 
of varying length. Considering the physical sig- 
nificance of the mentioned examples in multipar- 
ticle dynamics, it will be interesting to see how 
the proposed scaling behaviors show up in the 
experiments of the Third Millennium. 
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